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ANALYTICAL METHOD OF APPROXIMATING THE MOTION OF A SPINNING 
VEHICLE WITH VARIABLE MASS AND INERTIA PROPERTIES ACTED 
UPON BY SEVERAL DISTURBING PARAMETERS 

By James J. Buglia, George K. Young, Jesse D. 

Timmons, and Helen S. Brink worth 


SUMMARY 

An analytical method is presented which approxi- 
mates the flight-path deviation of a variable-mass, 
variable-inertia, spin-stabili zed vehicle under the 
influence of (a) initial pitching motions, (6) thrust 
m isalinenients, and (c) principal-axis misalinemeni , 
or dynamic unbalance . This method includes a 

first-order analysis of the effect of jet damping on 
the result in g flight path . 

An extensive analysis was carried out on an IBM 
70 j electronic data processing machine , by integrat- 
ing numerically the six-degree-of freedom equations 
of motion of a thrusting body with variable mass and 
inertia. This computer analysis was intended to act 
as a check on the method developed herein, and. also 
to present other investigators with exact numerical 
answers with which to compare their own results. 
The results from the analytical solution, are com- 
pared with the results from the numerical integration, 
and the results are found to agree well in all cases 
where no large change in inertia ratio was involved. 

The results are presented in terms of the angular 
deviation in space of the velocity vector from its 
undisturbed orientation, as a function of spin rate , 
for each of the disturbing parameters. It is believed 
that this velocity-vector deviation as a function of 
spin rate is a much more intuitive indication, of the 
dispersion than the body attitude in space, which is 
the most frequently used parameter in spin-stub dila- 
tion discussions. 

INTRODUCTION 

With the fairly recent advent of very-high- 
altitude rocket vehicles, the engineer has become 
faced with the difficult problem of adequately 


stabilizing his vehicle in the region where the 
atmosphere is so thin that conventional aerody- 
namic controls are rendered useless. One possible 
means of stabilization is the use of reaction con- 
trols in place of the aerodynamic controls. How- 
ever, in view of the diminutive size of some of the 
present-day probes and satellite injection stages, 
a control system of this type would impose a 
severe weight penalty on this final stage. Hence, 
some other means of pitch and yaw attitude 
stabilization must be considered. 

Ballislieians have long known that a spinning 
shell deviates from a predetermined flight path 
much less than a nonspinning shell. It thus 
seems natural to attempt to attitude stabilize a 
rocket vehicle by spinning it and, indeed, this 
method has beam successfully used already in the 
Vanguard, Kxplorer, and Pioneer vehicles, and on 
Scout and other vehicles currently being designed. 
Spin, however, is not a cure-all: spinning a vehicle 
gives it a constant attitude in inertial space, but 
not stability along the constantly changing tan- 
gent to tiie (light path. Hence, if a large change 
in flight-path angle is involved, the method of 
spin stabilization may not be applicable. Spin 
stabilization is very effective, however, in the case 
of a satellite launcher where the last stage, the 
injection stage, is final at nearly horizontal flight - 
path and attitude angles and it is desired to burn 
out this stage in such a manner that the final 
velocity vector will be nearly horizontal. 

Some investigators in this field (for example, 
refs. 1 and 2) have chosen as a criterion for spin 
stabilization the attitude of tin* body, defined by 
the two Killer angles in pitch and yaw. It is fell, 
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however, that a knowledge of the vehicle attitude 
alone is not the best criterion for determining the 
effects of various extraneous forces, because the 
attitude alone does not indicate the final orien- 
tation of the velocity vector at stage burnout. 
This velocity-vector orientation is believed to be 
a better indication of the dispersion effects sought. 
Consequently, the data presented herein are given 
as deviation of the velocity vector from its undis- 
turbed orientation as a function of spin rate for 
several external disturbing forces. 

The two methods in most general use for deter- 
mining dispersion ot a spinning vehicle, with the 
hoi ly attitude as a criterion, are the methods of 
Nicola ides (ref. 1) and Jarmolow (ref. 2). Xico- 
laides determines the dispersion in terms ot angle 
of attack a and angle of sideslip 0 his parameter 
is essentially % a 2 +>\ He also includes aerody- 
namic effects, which are neglected here. However, 
he considers the vehicle mass and moments of 
inertia, as well as the applied moment and foiwaid 
velocity, to he constants. His results, therefore, 
cannot be used to study the dispersion of a spin- 
ning rocket during its thrusting phase. 

damiolow, on the other band, considers the case 
of variable mass, moments ot inertia, and applied 
moment ; but he discusses only the body attitude. 
Also, analytical results are given only for the pitch 
and vaw rates in a body-axis system. Numeiical 
integration must still he used to find the pitch and 
yaw Euler angles. 

The method presented herein utilizes exponen- 
tial approximations for tfie variation ot the ratio 
of the applied moment to the pitch moment of 
inertia, and for tin* ratio of thrust to linear mo- 
mentum, which appears as a parameter in the 
translation equations. These, approximations 
make possible analytical solutions for pitch and 
yaw Euler angles, and also expressions defining the 
angular deviations of the velocity vector. I hese 
expressions include linear damping in pitch and 
yaw. The following additional assumptions are 
made: the ratio of roll inertia to pitch inertia, is 
constant: the vehicle is rigid and has rotational 
svmmetrv; then 1 is no roll moment; only small 
angles are considered; the motion takes place in a 
vacuum. This study was initiated to determine 
the effect of various disturbing parameters on the 
(light path of a solid-propellant rocket motor 
similar to the final stage of the Scout. 


SYMBOLS 


a, n, /),/•: 

constants 

<■' 

pitch and yaw damping coefficient, 
ft -lb-sec 

c= j> 1/sec 


constants 

H„ 

angular-momentum vector of ve- 
hicle with respect to center of 
gravity in /*,?/*, z b coordinate 
syst cm, slug-ft ‘feec 

H, 

angular-momentum vector of ve- 
hicle' with respect to center of 
gravity in x n y Cf z e coordinate sys- 
tem, slug-ftySec 

H, 

angular-momentum vector of ve- 
hicle and jet with respect to 
(‘(‘liter of gravity in x n y n z r coor- 
dinate system, slug-ft 2 /»ec 

i - j ^ ‘i 

unit vectors along the AY, 5 tr, 
and Z,-axis, respectively 

^ r . J i • ^ i' 

unit vectors along the A r -, ) 
and Z r -axis, respectively 

/v,/.-,'* 

vehicle moments of inertia about 
AY, )Y, and /,-axis, respec- 
tively, slug-ft- 

I 

J y or I x (where Iy Iyj, slug-fl- 

I y/.t Ay v 

vehiclc products of inertia about 
AY, V ft -, and //-axis, respec- 
tively, slug- ft" 

A, A,, A, A 

constants 

l 

length from vehicle center of grav- 
ity to nozzle exit plane, it 

M 

disturbing-moment vector, d/ Y i c 
+ Afrh-\- A/^k r , ft-lh 

M 

magnitude of disturbing moment, 
ft-lh 

m 

mass of vehicle, slugs 

HI H 

element of mass, slugs 

\t 

. lit. t 

U, 1/sec- 

/v r 

n—-j null! 

ans/scc 

]W 

angular-velocity components ol 
vehicle about AY, 1Y, and Z t r 
axis, respectively, radians /see 
unless otherwise indicated 

Te 

distance from center of gravity to 
jet exit, r ( A ir A ,vjrT .zk r , it 

r e 

magnit ude of r, , ft 
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t 

time, sec 

7 

t hrust , lb 

V 

vehicle velocity vector, ft /sec 

r 

magnitude of V, ft /sec 

V 

jot exhaust velocity, k^.\i r + 
V.'Y)c+V<‘.zK, ft /sec 

r. 

magnitude of V,, ft/sec 

j h 1 ,7 f) , - b 

coordinates measured along body 
AV, )V, and Z A -axis, respec- 
t ively 

A /, , l b 

body-fixed system of orthogonal 
axes, origin at body center of 
gravity 

j r 1 1fey-c 

coordinates measured along AV, 
>V, and Z,.-axis, respect ively 

A r , ) fJ Z e 

system of orthogonal axes which 
pitches and yaws with vehicle 
but does not roll, origin at body 
cent(*r of gravity 

♦( u >?/ u 

coordinates of element of mass 
causing dynamic unbalance 

CL 

angle of attack, radians 

it 

angle of sideslip, radians 

7 

orientation of velocity vector in 
vertical plain*, radians 

8 

thrust misalinement angle 1 , radians 
(except as noted) 

€ 

principal-axis misalinement angle, 
radians (except as noted) 

X 

orientation of velocity vector in 
horizontal plain*, radians 

(7 

total orientation of velocity vector, 
radians (except as noted) 

e 

pitch Euler angle, radians (except 
as noted ) 


yaw Euler angle*, radians (except- 
as noted) 

Oi c 

angular velocity of nonrolling sys- 
tem in inertial space*, a> }V jV- 
<jd z Mn radians/sec 

W}V,Wz«- 

Subscripts : 

components of w r along )V and 
Z,.-axis, respectively, radians/sec 

0 

denotes quantity at zero time 

X,Y,Z 

components parallel to AV, 3 V, 
and Z, ,-axis, re*spectivi*ly 

AV, Yv t Zc 

components parallel to AV, J V, 
and Z f -axis, respect ively 

A dot over 

a quantity denotes differential ion 

with respect to 

t ime. 


ANALYSIS 

The mathematical analysis of a spinning rocket 
differs in several respects from that of a spinning 
projectile. Jn the first place, the rocket lias mass 
and inertia properties which vary with time, 
whereas those of a spinning projectile are constant. 
Secondly, the rocket can have a moment due to a 
thrust misalinement. This moment, which is 
fixed to the spinning vehicle, gives rise to motions 
not found in the analysis of projectile motion. 
The motions due to an initial pitch or yaw rate 
and a dynamic unbalance an* common to both 
types of bodies, but t lie analysis becomes more 
involved when a variable-property vehicle is dealt 
with. 

The (‘({nations of motion programed on the 
IBM 704 electronic data processing machine for 
comparison purposes in t his st udy are t lie standard 
six-degree-of-fieedom (‘({nations found in most 
aeronautical textbooks (for example, ref. d ). 
The ('({nations programed, along with tin* various 
transformations that wore used, are presented for 
reference in the appendix. 

An approximate analytical method of predicting 
the deviation of the velocity vector due to the 
effects of thrust misalinement, initial pitch rate, 
and principal-axis misalinement will be developed 
in this section. Two coordinate systems will be 
used in this analysis. The first is a 

body-fixed coordinate system with axes fixed at 
the center of gravity of the vehicle and the 
coordinate x ft extending along the body longitudinal 
axis. The second coordinate system (x ef y e ,z f ), in 
which the equations of motion will lx* derived, is 
referred to a body-fixed system of axes which can 
pitch and yaw with the vehicle but does not roll 
wit h it . The t wo sets of coordinate axes are alined 
at /. — (). (See fig. 1 .) The pitch Euler angle 0 
and tlie yaw Euler angle \p are assumed to lx* small 
so that the older of the Euler rotation is immate- 
rial. This assumption of small angles simplifies 
the integration for 7 and X, the flight-path angles 
in the vertical and horizontal planes, respectively, 
and also permits these flight-path angles to 
combine veetorially to give <r, the angle between 
the velocity vector and some fixed reference line, 
taken here as the longitudinal axis of the vehicle 
at zero time. 

The angular-momentum vector of the vehicle 
in the system instantaneously fixed in 




VACUUM MOTION OF SPINNING VEHICLE WITH VARIABLE MASS AND INERTIA PROPERTIES 


O 


of gravity is 



(:*) 


The rate of change of angular momentum with 
respect to t ho vehicle center of gravity including 
jet effect is 


r/H r <}\{ r , v __ , . . WT v 

~~/j~ : : ^ ■ + w f xH f + w r r X (w, X r f — \ ,.) (4 ) 

where w f XH ( . takes into account the rotation 
of the coordinate system. 

It is usually conventional in missile work to 
choose the body axes as principal axes and to 
assume a symmetrical vehicle, so that 



With these assumptions, equation (1) may be 
written in the nonrolling system as 

H f l x pl, + /w,vj f d -IwzMc (t>) 

If the further assumptions are made that the 
jet is located on the AVaxis and the jet exit 
velocity relative to the nozzle exit is directed 
along the negative AVaxis, then 



Expansion, of equation (A), with the aid of sub- 
stitutions from equations (2), (4), and (0), yields 

( / Y V "f" f Xp) if ( Iw V' ( V /w }' P + / K ]) (&z r. “T M r i « Yc ) j r 

(Ic o Xc ~r Ia Xe — I x p(^y c ~\~ th co^ f )k c -- M (9) 

If damping in roll is neglected (that is, if f x -()), 
e([uation (9) becomes, in component form, 

I x p—J\ f x (10a) 

l^Yc~^ I x )^zr J rC' Oi Yc M y (10b) 

l^Zc I X P^y c~T C &Zc.— ( 10c) 


where 

c f = mr 2 e J ci 

The quantity c' jl(—c) will be considered constant 
for this analysis. 

If the assumptions are made that the disturbing 
moment lies in the vertical plain* at zero time 
and that there is no moment about tin* roll axis, 
1 hen 



Vv = 0 'J 

My— M cos pt V 

M x —M sin pt J 

(id 

and (‘({nations (10) become 



l x p-i) 

(12a) 

Ja Ye 

X /v/^f + c'Wy c - M (“OS pt 

(12b) 

I*ZC 

— I xpay c -\- C r a Xc ~ \l sin pt 

(12c) 

E({uation (12a 

) immediately shows that 



p^ ( \mstant j> 0 


Let 

lx 

~j~ Po~ ft 

(l:;) 

and 


(14) 


It is assumed that the ratio l x /I is constant; 4 f 0 
and k an* also constants, found by fitting an 
exponential to tin* function Mil as shown in 
(‘({nation (14). Then equations (12b) and (12c) 
become 

COy c -\- n<ti%c-\-Ca)y c —- M 0 e kt cos pj (I mi) 
uzc~ n-a Yc -\-cw Xc — M 0 ( jkl sin p (i t (lob) 


Thes(* (‘({nations can be integrated simply, and 
with the initial conditions «jv(0) “ aqv.o, a Zf (())= 0 
give 


( \)e ct cos ni— ( \v a sin ni 
cos pJ~r (']('** sin pj 

a%c— (uyc. <j — ( Vic' 1 ’* si 11 ni- (\t~ ct cos lit 

-- ( V' sin pj — ( \e kt cos pj 


where 


, M t (p„ — n) 

1 (C-\-ky 2J r- (p„— /(,)" 


( 1 On) 


(10b) 


(17a) 
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.\fjc±k) 

(c + k ) 2 + (p„ — n)~ 


(17b) 


Now, bv virtue of the original assumption that 
the ])it(‘h and yaw motions are small, and as a 
consequence of the coordinate system used, the 
following expressions can be written: 

d—J'uy c f it (lSa) 

i p— S (isb) 


Substitution of equations (16a) and (1Gb) into 
these integrals with the initial conditions 
0(O) = ^(O)™O gives 


e—A(e kt cos pJ—\)-\-Iie kt sin pj 

— I){e~ ct (‘os nf—\)-\-Fe~ ct sin nt 

\p~ A( kt sin p 0 t—H(e kt cos p„t — 1) 

— I>(~ ct sin ni — FA(~ ct cos nt — 1) 

where 


.1 = 


IF 


_ < ik- < \p, 

k 2 f pi 

( 'iPo- X r ( \ k 

fr \ pl ' 


I h 


Cia'iv.,, — f V) — n('i 


(19a) 

(19b) 

(20a) 

(20b) 

(20c) 


K 


— ( A 1 c( 
c 2 -\ n 1 


(20d) 


Those equations give time histories of the body 
attitude angles in an inertial coordinate system 
taken to be coincident with the body coordinate 
system at zero time. These angles 6 and p are 
the angles commonly used to discuss the effects of 
spin on the stability or dispersion of a vehicle. 
However, body attitude alone does not adequately 
describe the complete motion in space. The act ual 
dispersion of the vehicle, defined herein as the 
angle between the velocity vector at zero time and 
the velocity vector at burnout, is a much more 
meaningful parameter because it depends on the 
translation of tin 1 vehicle, and hence on the quan- 
tities which define its translation history, as well 
as on the rotational parameters. The preceding 
argument can be substantiated from physical 


considerations. The equations defining the rota- 
tional Ji iot ion of the vehicle are completely inde- 
pendent of the translational motion. (See, for 
example, eqs. (Al) to (AG).) Hence, for given 
inertia characteristics, spin rate, and disturbance, 
the rotational motion of the vehicle is completely 
defined. Now, in order to get the actual transla- 
tion in space, the body attitude parameters must 
be combined with the translation parameters; 
for example, equations relating quantities such as 
thrust and mass ratio are integrated to give the 
translational motion. As an illustration, it is 
certainly possible to have a spinning vehicle point- 
ing vertically upward and moving horizontally at 
ignition. Then as the motor fires the velocity 
vector climbs above the horizontal. Thus, the 
final direction of the velocity vector has (‘banged — 
the dispersion as defined here is reduced — but the 
body attitude has remained constant. An effect 
such as this one is not evident when only the rota- 
t ional motion is defined. 

In order to determine 7, the orientation of the 
velocity vector in the vertical plane, reference is 
made to tin* normal equation of motion of figure 
2(a). Since a — 7, if only small angles are 

taken ii to account this equation becomes 


<h , T T 
,/rm V y mV 



cos 


pj 


( 21 ) 


In a sii filar manner a differential equation for X, 



(a) 


mVy - sin 8 cos p Q t cos a +T cos 8 sin a 
mVy = 78 cos p 0 t + 7{8-r) 

(a) Vertical plane. 

Fho kk 2. Filiations of motion. 
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the* orientation of the velocity vector in the hori- 
zontal plane', can lx* obtained from figure 2(b): 


</\ , T x T , , T . , 


( 22 ) 


Since only the changes in 7 dues to external dis- 
turbances rather than the natural change in flight 
path due to gravity are* of interest, the* (*ffeot of 
gravity is omit te*d. 

Now, in order to integrate equations (21) and 
(22), the* quantity TjmV must be put into a form 
that is easily integrated. It can be shown that a 
plot of Tim l ” as a function of time can be* approxi- 
mated to a suflieient elegree of ae*ouiae*y by an 
exponential function. Thus, define 


T 

m\ 




(24) 


an el 



T cos 


/- h 

k. 


In the pre>cess of integrating (‘((nations (21) and 
(22), the quantity appears in several in- 

stunces. Therefore, to further simplify the in t e»- 
m-ation, it was assumed that 






(24) 


mV \ - T sin 8 sin/^/cos £+ T cos 8 sin 
m V X = T 8 sin p Q t + T $ 

(h) Horizontal plane. 

ITeaJKK 2. Concluded. 

If equation (2.4) is substituted into equations 
(21) and (22) along with t lie expressions for 6 and 
\p from equations (19), integration of the resulting 
linear different ial equations, utilizing eqnut ion (24) 


and the initial conditions 7 ( 0 ) A(0) 0, yielels 


7-c 


T( ^ k p f r Ae ^ kl ' (( \+lr) -p„ Be '<■’* » k] ' he v * 1 ( ^ Cl p 0 e ■'*+ k) ' + Be ' k} 1 (( -(•• k) 

1 1 lL ((?>-\kY+pt + n+pd <I()S/ + L {(~ 2 +ky+pi 

t 5 eWpol . . r/)(( \-c)e {r '*- eU + Km®*- e 'n rnDe^- cU -Ke^~ c)l {r,-c)~\ . 

Ci+nZJ L w 3 - c )*+ n * J L (a-cY+n* J 


sm nt 


~Af« Df'-n V A {V t +lc)-B Vv IHCy-c) +Eu_A D <SO, ~l\ , 

- ('■: (\ J L (Ci+W+Pl ((\-c)*+n* (Y (Y WV 


(CU+ky+pl {(’i-c)*+n' 

x= <>-'*** kc P M e « i ( . os 


p„Be tc - +k) 

(rvn-fYpi 


- r .s 

riky-vri «'1+/»rJ L {(\ • kf 

, ■ t TnIh , ~*- rt '-(T\,-c)Ee'~2- r >n ri)e l ~r fU (F,- c) + /i/',Y ,r V '"1 

+_- sin i> 0 t+\ I cos nt — I - _* s 

<1-1 piS L (c.-r)* | »* J L J 


sin nt 


r-ivyfr’,-) k) _ j,j> nv- (< h-c)E, k /n\ (2 

L c, + (’■. J L (/ 2 i ky+pi r - 1 C {('.,-cyy «- rj/ 


(26) 


023963 -02 
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Til pm, because of the assumption of small angles, 
7 and X can be combined vectorially to give a, tin* 
angle between the fixed reference line and the 
velocity vector: 

o- - vV+A 2 (27) 

A disturbing parameter consisting of a princi- 
pal-axis misalinement can be bandied in tlie fol- 
lowing manner. Assume that the mass clement 
causing tin* dynamic unbalance' can be placed oil 
the center line of the vehicle a distance \j u from 
the center of rotation, and that it has a mass 
m u , as in sketch 1 . The moment about the center 



Sketch 1. 


of gravity due to this mass element is 

\f pU' u y H m u (2S) 

Now, tlu* product of inertia due to this mass is 

I xv -■-*«?/ (29) 

Hence, the moment can l>e written in terms of the 
product of inert ia: 

M~pi1 xy (:«)) 

It may be desirable to express the product of 
inertia in terms of the principal-axis misaline- 
ment angle e. From the properties of moments 
of inertia it is known that the angle e between 
the principal axis and the axis about which the 
moments and products of inertia are computed 
is given by 

lan2e-= : Ixr . Oil) 

* X 1 

If, as will be the usual case, the product -of-inert ia 
term is due to construction technique or placing 
of instrumentation, the product of inertia will be 
constant, the small change due to the oenler-of- 
gravity shift during burning being neglected. In 
this case, e will be a function of time. Tin* time 
history of e can bo computed from equation (:U) 
since / v and I are known at any time. Then 


with the assumption c<< 1, substitution of equa- 
tion (HI) into (‘([nation (HO) gives 

'/ -/«(/)( 1 -/) (:«) 

If this expression is used in equation (14), and \f„ 
and k art* found by a curve (it , solution of equa- 
tions (25) and (2b) will yield histories of 7 and X 
due to a principal-axis misalinement e. 

It should be noted that, to the degree of approxi- 
mation used herein, the (‘fleets of the various 
disturbing parameters can be combined linearly 
to yield the total dispersion. For example, equa- 
tions (25) and (2b) may lx* written as a linear 
combination of terms each of which is a function 
of o>yc,i , d/„, or 8 but not of any two of the variables. 

It should be further noted that, again to the 
degree of approximation assumed here, the angles 
7 and X, and consequently the total dispersion <r, 
vary linearly with o) Yco and A/„, if taken sepa- 
rately, so that it is not necessary to make a. 
separate calculation for each value of moment or 
pitch nite desired. This linearity is supported 
by tin calculations performed on the IBM 704 
computer, as illustrated in several of the figure's. 

RESULTS AND DISCUSSION 

A comparison of the results of the numerical 
integration (on the IBM 704 computer) of the 
exact '‘(pnit ions of motion given in t lie appendix 
with the results of the analytical solutions de- 
veloped herein is made in figures 4 to 1 b. The 
quantifies 6 , 7 , X, and a were calculated from 

equations (1 9a), (l Ob), (25), (2b), and (27), 
res peci ively. The constants used arc given in 
tables 1 and II. Figures 7 to lb are plotted for 
/ = BO seconds (burnout). These figures serve a 
double purpose. Since not everyone has access 
to an IBM 704 or to a computer of sufficient 
capability and storage capacity to permit the 
progra ning of a problem as complex as a six- 
degree -of-freedom rigid-body analysis, many in- 
vest ig» tors have probably developed analytical 
results without being able to check the accuracy of 
their methods. Thus, the figure's presented can 
serve to check these methods and indicate regions 
of validity and applicability. At the same time, 
the figure's serve to illustrate tin* accuracy of the 
analytical method developed herein. As can he 
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TABLE I. PHYSICAL AND PERFORMANCE 
CHARACTERISTICS OK ROCKET MOTOR 
ESEI) IN CALCULATIONS 


Roll moiiK'iit of inertia, si ug-ft 2 .. 5.2o 0.1^ 

Pitch (yaw) moment of inertia, slug-ft 2 .. 84.15 0.88/ 

Mass, slugs _ ..... _ 18.05 0.18/ 

Applied moment , ft -ID (5 in radians) 10,7105 : 1505/ 

Thrust, lb ... .. . 8,000 

Burning t ime, sec . . 80 

Initial velocity, ft/aec _ . . 10,000 

Final velocit y', ft /sec _ 10,810 

Computed constants (for 0. 10° thrust misalinement ) : 

0.01588 

k>. . . . . 0.0200 

I _ - 0.008 

(\ ....... 1.8 80 

r, . 0.050 

5, radians 0.001745 

c oy e . radians/sec . .. , . ...... 0 

/ v// - -- 0.104 

.U„. 1/sec 2 ... . 0.5470 

k ... _ 0.0502 


TABLE II. PHYSICAL AND PERFORMANCE 
Cl I A R.ACT ERISTICS OF SPHERICAL ROCKET 


MOTOR 

Roll moment- of inert ia, slug-ft- 888.0 10.0/ 

Pilch (yaw) inonienl of inert ia, slug-ft 2 408.0 10.0/ 

Mass, slugs 10.88 0.51/ 

Thrust, 11) 8.400 

Burning I ime, sec . _ 80 

Moment, ft-lb 5 | 0.01/ 

Initial velocity, ft, sec 10,000 

Final velocity, ft/sec 20,050 

(Computed constants: 

k i 0.02125 

k 2 . .. . 0.0700 

Av _ 0.200 

7\ 1.280 

C 2 - .... 0.107 

5, radians . .. . . _ . 0 

Wyc.o, radiuns/see ...... 0 

i x :i . ... __ o.o8 

M„, 1/sec 2 . . . _ 0.0124 

k ....... .... . 0.042 

soon from a general perusal of figures 3 to 15, t lie 


accuracy obtained is quite good in the range 
invest igated. 

Figure 3 shows the pitch Euler angle 0 [dotted 
against, the yaw Euler angle \f/, for the last few 
seconds of burning of the rocket motor, under the 
action of a moment produced by a thrust misaline- 
mont 5 of 0.10°. This figure shows, essentially, 
the attitude of the body in space during this time. 
It can be set'll that, although the analytical method 
predicts frequencies of nutation and precession 
which are slightly high, the amplitudes of both 


0,deg 



Fit; no; 8. Undamped body attitude. }>„ 80 rpm; 

5-0.10°. 

motions, which are the parameters of primary 
interest, are predicted rather well. This curve 
dot's not include the effect of jot damping. 

Figure 4 shows lime histories of the deviation 
of tin 4 velocity vector a. The analytical method 
predicts results which are slightly conservative. 

In figures 5 and b are plotted time histories of 
the pitch Euler angle d and the yaw Euler angle \j/ y 
respectively. Here it will lx* noted that, while 



0 4 8 12 16 20 24 28 


Time, sec 

Fkukk 4. -Deviation of velocity vector as a function of 
time, for a thrust misalinement angle 5 of 0.10°. Un- 
damped body; p a — 80 rpm. 















APPENDIX 

EQUATIONS OK MOTION IN BODY-FIXED AXIS SYSTEM 


The symbols 
as follows: 

used in this appendix are defined 

FxJ'yS'z 

components of external force along 
the body A-, T-, and Z-axis, 
respectively, lb 

7 A'./ yJz 

moments of inertia about body 
A*-, Y-, and Z-axis, respectively, 
slug-ft 1 ’ 

7 xz 

product of inertia, slug-ft 2 

I 

I Y or l z (where 1 Y —I X ) } slug-ft 2 


direction cosines betwe(*n inertial 
axes and body-fixed axes, 7=1, 
2, or )> 

M x JfyM Z 

external moments about body AT-, 
Y- % and Z-axis, respectively, 
ft -lb 

in 

mass of vehicle, slugs 

7h77* 

angular-velocity components of 
vehicle about body A”-, K-, and 
Z-axis, res])ectivelv, radians/sec 

r n X 

longitudinal distance from vehicle 
center of gravity to jet exit, ft 

inrmu 

t rans format ion mat rices 

u,rjr 

velocity components of vehicle 
along A-, T-, and Z-axis, rospoc- 
tively, ft /sec 

V 

vehicle velocity vector, ft /sec 

XX,Z 

orthogonal body-axis system 

yj/x 

coordinates of center of gravity of 
vehicle along A-, } and Z-axis, 
respectively, ft 

xx,z 

fixed orthogonal axis system, coin- 
cident with body-axis system at 
zero t ime 

7 

orientation of velocity vector in 
vertical plane, radians 

6 

pitch Eulor angle, radians ((except 
as noted) 

X 

orientation of velocity vector in 
horizontal plane, radians 

<T 

total orientation of velocity vector, 
radians (except as noted) 


(j) roll Eulor angle, radians" (except as 

noted) 

\j/ yaw Euler angle, radians (except 

as noted) 

A dot over a quantity denotes differentiation 
with respect to time. 

The rigid-body equations of motion as pro- 
gramed on the IBM 704 electronic computer are 
the standard six-degree-of-freedom equations of 
motion found in most aeronautical textbooks on 
dynamic stability (for example, ref. eqs. (10-1 1)), 
with damping in pitch and yaw added to the 
second and third moment equations, respectively. 
These (‘({nations are repeated here for reference. 
The tin ee force (‘({nations an* 

m ( it — rr+w<j) — X Kx (A 1 ) 

m ( v -- wp + nr) — X7 y V ( A 2) 

m (ib— u ({ + rp) — X Vz (Ad) 

and the moment equations including jet damping 
are 

I. p—Ix/j \- ( Iz — I y) <//' — IxxP<i~ X^ Ay ( A4i 

7 y(f - j (7 Iz)rp-r IxAp 2 — !*) 

-f - 1 y <i + mr; t x q — X d / y (A 5 ) 

lyj'— /v 4- < / r — /.V ) l><l+Ixz<P' 

1 w/t. y /* yih (AOj 

r Fhe [uantities used hero are referred to a set of 
axes ib ed to and rolling with the body, whereas 
the analysis in the body of this paper is referred 
to a set of non rolling axes. 

The quantities u, i\ w and p, (/, r are found by 
integration of these (‘({nations by tin* Runge- 
Kutta net hod, and are then transformed to a set 
of iner ial axes. 'Die following assumptions are 
made t ) simplify tin* integration: 

M >ments and products of inertia vary 
in early with t inn* 


IS 
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Muss varies linearly with time 

Damping coefficients (mrj tX ~\-I) are constant 

The inertial quantities are found by means of 
the following transformation equations (ref. 4): 


where the transformation matrices [1 \ ], [7 2 ], and 
are given hv 


[TA-- 


i U 


cos \p — sin t p 0 
sin \p cos \p 0 

0 0 I 


cos 6 0 
0 1 
-sin 6 0 


sill 6 


(AS) 


( AO) 


cos 8_ 

't 0 0 

1 7 :{ ] — | 0 cos 0 —sin </> | (Alt)) 
_{) sin 0 cos 0 J 

Multiplying the matric(*s gives (‘<juation ( A 7 ) as 

/j m i u 

If | | /_> Hi 2 ih i' I (All) 

_/ :i /oj L ?r . 


whei 


/ j — cos yp cos 0 
Ap — sin \p cos 6 
/ ;i — —sin 8 

///, — cos yp sin 6 sin 0 — sin \p cos 0 

/// 2 — sin \p sin 6 sin 0-j cos yp cos 0 

///*= cos 0 sin 0 

ti j = cos i// sin 0 (*os 0-( sin yp sin 0 

/o = sin \p sin 0 cos 0 — (*os xp sin 0 

//{—COS 0 COS 0 


> (A 12) 


The angular rates y/, (/, /* about the body axes 
are transformed to angular rates 0, 6 } \p about the 
space axes by means of the following transforma- 
tions (ref. 4): 


(A 13) 


Application of these transformations to the body- 
axis quantities u } i\ v s and y/, (/, /’ yields time 
histories of 7, 7/, ™ and 0, 0, yp. Then, from sket(*h 
2, 7 , X, and <7 can be found: 


.7 


M 


0=y/-|-tan 0(<y sin 0 + /* cos 0) 

// 

=r/',i[7’ 2 j[7’,i 


(A7) 

0 = -(j cos 0 /* sin 0 l 

0—- --- (« sin 0 + /* (‘os 0) 
cos 0 

- r - 






Sketch 2. 


Thus, 


7 — tan 1 , ■ <j ' 

X— tail -1 ■! / 


(A 14) 


cr — t a 1 1 


-i d-.//‘ 
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